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Abstract 

We present a self-consistent solution of the finite temperature gap-equation 
for A<3? 4 theory beyond the Hartree-Fock approximation using a composite 
operator effective action. We find that in a spontaneously broken theory not 
only the so-called daisy and superdaisy graphs contribute to the resummed 
mass, but also resummed non-local diagrams are of the same order, thus 
altering the effective mass for small values of the latter. 
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INTRODUCTION 



The study of $ 4 - theory at finite temperature is of great interest for a wide field of 
applications. A self-interacting scalar field serving as a simple model for the Higgs particle 
in the standard model of electroweak interactions may allow the study of symmetry changing 
phase transitions. In fact, the order of the electroweak phase transition plays a crucial role 
in the framework of cosmological scenarios as well as for the badly understood process of 
baryogenesis [|I|]. Despite the simplicity of the model, one may at least hope to gain some 
insight in the mechanism of the phase transition. Moreover, the theory is a suitable test 
ground for analytic non-perturbative methods, e.g. variational methods [|| as well as for 
lattice simulations |J. 

High temperature symmetry restoration in a spontaneously broken theory was already 
noted by Kirzhnits and Linde Q and worked out quantitatively subsequently ||. The 
convenient tool to study the behavior of the theory turns out to be the effective potential. 
Whereas the critical temperature is already determined by the one loop potential and thus 
relatively simple to find, the order of the phase transition depends on the detailed shape of 
the potential which requires also the analysis of higher loop contributions, even for small 
coupling constant. 

In particular, one finds, quantizing the theory around the classical non-trivial minimum 
that the one loop self-energy at high temperature behaves like (mr) 2 ~ {-\f\T) 2 , indepen- 
dent of external momenta and the value of the classical minimum. Thus, for large T the 
thermal mass dominates over the tree-level mass and the minimum of the effective action 
becomes the trivial one. The theory exhibits two important features: non-temperature sta- 
ble vacuum and effective temperature dependent mass. It was consequently proposed 
to move the tree level mass into the interaction part and to start perturbation theory with 
the free Lagrangian + fi 2 )$ including a yet undetermined mass parameter Q. Then 

one calculates the effective potential and fixes the parameter by the 'principle of minimal 
sensitivity' dV{$ min )/dVL = 0. 
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A more systematic approach of a self-consistent loop expansion was developed already 
some time ago by Cornwall, Jackiw, and Tomboulis (CJT) in their effective action formalism 
for composite operators ||. The basic idea is to introduce a bilocal mass operator in the 
generating functional instead of the local mass Q. Then one defines the generalized effective 
action T as the double Legendre transform of the generating functional which is now not 
only a functional of the expectation value of the field but also depends on the expectation 
value of the time-ordered product T$(x) $(?/). The principle of minimal sensitivity gets 
replaced by the the so-called gap-equation 5T(^,G)/5G(x,y) = which is employed to 
eliminate the exact two-point Green function G(x,y). 

This paper is dedicated to investigate on the solution of the gap-equation beyond the 
Hartree-Fock approximation. Whereas the latter amounts to simply include only a local 
mass term in the ansatz for the Green function |7||§, thus being a slightly more elegant 
reformulation of the variational approaches mentioned above, we will consistently || include 
also the non-local contribution appearing in a spontaneously broken theory. 

Our approach is as follows. We split the inverse Green function into a local mass term 
absorbing the usually resummed graphs, and a non-local self-energy part. By this ansatz, 
the gap-equation is recast into a non-linear integral equation for the self-energy, which we 
solve approximately by a suitably chosen ansatz. Self consistency requires the unknown mass 
parameter to satisfy an 'effective-mass' equation, the latter differing from the one found for 
the effective mass in the superdaisy resummation by an additional term which dominates in 
the limit of vanishing mass parameter. We discuss the behavior of the mass parameter near 
the critical temperature. We lay great emphasis on dealing with the divergencies correctly. 

CJT COMPOSITE OPERATOR FORMALISM 

It is useful to briefly review the basic steps in the construction of the generalized effective 
action for a classical action /($) 0. One starts with the generating functional in Euclidean 
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space-time 



Z(J,K) = J D©exp- J(©) + J J r (§)©(§) + 



OG 



©(§)£(§, t)©(t)) 



§ t 



and W(J, K) = log Z(J, K). After defining the classical field <f>(x) and the Green function 
G(x,y) by 

5W(J,K) 5W(J,K) 1 



(x) = <$(*)) 



■ t (G(x,y) + <t>(x)(t>(y)) 



1) 



SJ(x) " rv " v ^^"J,*> 5K(x,y) 2' 

one performs a Legendre transformation to find the generalized effective action 

T(cf>, G) = W(J, K) + J (j>(x)J{x) + l -J J (<f>{x)K(x, y)<f>(y) + G(x, y)K(x, y)) 

x x y 

where J,K at the r.h.s. have to be taken at the solution of Eq. ([!]). From the above 
definition, one immediately deduces the inverse relations 



54{x) J(*) + jK(*>vMy)> 5G{Xjy) 2 

The last equation tells us that if one wants back the standard effective action defined by 
K(x, y) = one has to impose the so-called 'gap-equation' 

sm g) _ 

5G(x,y) • 

In order to give a series expansion of T(4>, G) we define the functional operator 



ST(</>,G) 1 



K(x,y). 



8<j>(x)5<j>(y) 



and one finds 



r(0, G) = I(<P) + ^Tr HDoG- 1 ) + ^MD-'G - 1) + (0, G), 

D (x, y) being the free propagator derived from the part of the action which is quadratic in 
the fields. The quantity T^ 2 \(j),G) contains all two-loop contributions and higher and has 
to be calculated as follows: Shift the field $ in the classical action by 0. Then J($ + 0) 
contains terms cubic and higher in $ which define the vertices. r^(0, G) is given by all 





(a) (b) 
FIG. 1. Feynman graphs contributing to the lowest order self-energy. Lines denote the full 

Green function G(x,y), the tree-vertex is given by A</>/3! and the four-vertex by A/4! . 

two-particle-irreducible vacuum graphs with the given set of vertices and the propagators 
replaced by G(x,y). The gap-equation reads 

5T^((f), G) 



G-\x,y) = D-\x,y) + 2- 



SG(x,y) 



A$ 4 THEORY 



The self-energy graphs we are going to resum are shown in Fig. 1. Although the non-local 
contribution Fig. lb is formally of order A 2 it contributes to order A in the false vacuum since 
already the tree-level value of the expectation value of the field is ~ mj y/X. Consequently, 
the three- vertex A0/3! is of order a/A which means that in the case of a spontaneously broken 
theory the non-local contribution (b) is as important as the local one (a). We will show that 
this is also true at finite temperature. Thus, the Hartree-Fock approximation involving only 
the local graph which corresponds to the superdaisy resummation at finite temperature is 
incomplete. Including the non-local graph Fig. lb in the self-energy, the gap-equation in 
Fourier-space reads 

M 2 

, " > 

G~\k,u) = P + „ 2 + (-m 2 + ^) + \TY.j (03<?(a/, V) + 

+ ^) 2T E / -§^ G ^^ G ^ ~ % ( 2 ) 

where we use imaginary time formalism with u = 27inT, n G Z. The first sum/integral 
corresponding to graph (a) does not depend on external momenta and energies and thus can 
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be treated as mass term. We employ the ansatz 

G~ l {u), k) = k 2 + u 2 + fi 2 + U(u, k) (3) 

and subsequently absorb all local contributions in the mass parameter fi which has to be 
determined self-consistently afterwards. Plugging (Q) into @, the gap-equation transforms 
into an equation for H(u, k), 

d 3 p 1 



nM) = M 2 -,i 2 + i(A0) 2 T£ / 



(2tt) 3 p 2 + u' 2 + fi 2 + U(uj',p) 
1 



x 



• (4) 

(k - p ) 2 + /i 2 + (u - lu') 2 + TI(lu -u',k-p) 

Simple power counting now reveals that the expression involving the integration is logarith- 
mic divergent for large momenta p. To isolate these divergencies, we subtract and add a 
term of the form 

r d 3 p 1 1 f ,x 

yJ (2rr) 3 p 2 + uj' 2 p 2 + oj' 2 + /j? 1 ' 

which in the high temperature (T ^> jj) limit evaluates to 

' V, V= ^(l- 7 + log2vr + log^), (6) 



4vr/i ' 8tt 2V ' b to A 

where we introduced a UV-cutoff A. The difference between the last summand in Eq. (|J) 
and the subtraction (||) is a UV-finite quantity. Consequently, in the high temperature 
limit, we are allowed to pick only the zero-mode in the thermal sum in the difference (f|) - 
(H). The lo' = contribution of the subtraction exactly cancels the first term in (|6]), leaving 
us with |A 2 2 X> as divergent part of the r.h.s. of the gap-equation It is now possible to 
absorb the logarithmic divergence in the mass renormalization and we choose the parameter 
H such that 

\ 2 rh 2 

M 2 -fi 2 - ——V = (7) 

which renders the self-energy Tl(u), k) in equation (f|) a finite quantity. Furthermore, one may 
simplify equation (|J) by considering the to = mode of the self-energy only since the non- 
static modes are at least suppressed by a factor 1/T 2 . Thus putting n(o; = 0, \k\) = U(k) 
we end with the integral equation 
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oo 1 

n(jfc) — k f dpp 2 1 — (8) 

i i {p 2 + 1 + fl{p)){p 2 + k 2 + 2pkz + 1 + il{\jp 2 + k 2 + 2pkz)) 
where we introduced the dimensionless quantities p = /ip, U(p) = fi 2 fl(p) and the parameter 
« = T(A0) 2 /(87rV)- 

SOLUTION OF THE GAP EQUATION 

Without going in the details of solving (Rl) one can already make some useful statements 
about the behavior of fl(k). Firstly, the denominator of (|8|) is a positive definite quantity, 
thus ll(k) > 0. Furthermore, assuming that the solution is bound, i.e. fl(k) < fl ma x, the 
denominator may be estimated form below (above) by (p 2 + l){p 2 + k 2 + 2pkz + 1) and 
(p 2 + 1 + fl m ax)(p 2 + k 2 + 2pkz + 1 + flmax) respectively. Carrying out the integration on 
the r.h.s. in (||), one finds the following bounds for fl(k) 

7T k n(fc) 7T k 

— arctan — = < < — arctan — (9) 

k 2^Jl + fl max K k 2 

which in turn determine the asymptotic behavior, fl(k) = n 2 n/(2k) + (D(k~ 2 ). 

Secondly, we note that due to the z-integration, odd powers of z vanish in the integrand, 
which in particular means that tl'{0) = 0. 

Using this property and the asymptotic behavior, we are already able to approximate 
n(^) by a suitably chosen function. Due to the non-local character of the integral equation, 
however, we have to know the function Tl(k) for all momenta even if we were only interested 
in the low momentum behavior of the self-energy. Consequently, each kind of series ansatz 
is not a clever choice since it can never exhibit both, the correct small momentum behavior 
and the large k limit. We thus employ a Pade (rational function) approximation which 
will turn out to be a remarkably good global approximation for n. The simplest ansatz 
consistent with fl'(0) = and the asymptotic behavior reads 

A 2 + ^Kk A / \ 

n (k = T ^-, A = A(k), a = a{n) 10 

a 2 + §fnk + k 2 




10 20 30 40 50 ka PP a 

FIG. 2. Solution of the integral equation for II at k = 10. The dots represent values found by 

a 64-point Gaussian quadrature of the exact integral equation. The lower rigid line corresponds to 

the analytic approximation of the self-energy employed throughout this paper. We also indicate 

the exact asymptotic behavior of II (upper line). 

and a, A are functions of k which we have to determine such that the approximation fits as 
well as possible the true solution. 

From the form of the denominators in the integrand (|j) it is evident, that fl(p) contributes 
most to the integral for small values of p since for large values of p, IT being a decreasing 
function is neglectable compared with p 2 . For this reason we choose the ansatz such that 
(|) is fulfilled exactly for 11(0) and I1"(0). Plugging (|I"0"D into the integral equation, one 
could — at least in principle — carry out the integration which results in two complicated 
transcendental equations for a, A which have to be solved numerically anyway. We, however, 
will not do so, but merely discuss the limits k — > and k — > oo. It turns out (see appendix) 
that for large values of k, the coefficients exhibit a power law behavior, namely A(k) ~ 
Aq k 2 / 3 , cl(k) ~ ao ft 1 / 3 with the constants Aq = 4, 8611 . . . and ao = 4, 1140 . . .. For small 
k we find A(k) ~ \[§ixk v I 2 and a(n) ~ \^12 + . . .. Values between the small n and large k 
region can be found by numerical integration, see also the table given in the appendix. In 
Fig. 2 we compare the approximation ([To|) with the exact one found by numerical integration 
of Eq. @. 
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The solution for ft may further serve to eliminate the parameter M as defined in (^j), 
which amounts to evaluate the trace of the propagator, 

TYf^P- l - (11) 

where as before, we have to extract the UV-divergent part prior to be able to make the 
static approximation u = in the sum. Again subtracting and adding the UV-dominant 
contribution 

v r d 3 p 1 _ T 2 /jT A 2 2 

\ J (2vr) 3 p 2 + ^ 2 + /i 2 ~12~4^ + 8^ 2+/ " ' ( ' 

(0) " dH) is a UV-finite quantity, and we are thus allowed to drop all modes except for the 
zero-mode in the thermal sum. The difference may be expressed by the function 

" f{k) = * I V + /»')&»■ + <»' + n(p)) (13) 

which can be calculated analytically using the approximation for fl together with the known 
parameters a and A. We mention that, as before, also this function exhibits a power law 
behavior for large values of the parameter k, F(k) ~ F k 1 ^ 3 , F = 0.9963 . . ., and goes like 
F(k) ~ ktc/6 for small k. Finally, putting together the definition of M with the condition 
(0), one can transform the gap-equation (|2|) into an equation for the mass parameter /j,, 



/i 2 + -^-V = -m 2 + -0 2 + - ( — - ^- (1 + F(k)) + — + /i 2 P ) , (14) 



which is nevertheless meaningless until we absorb the divergent contributions A 2 , D in a 
redefinition of the bare quantities m, A and 4>. 



RENORMALIZATION 

The renormalization of A$ 4 -theory is delicate task ||10|| . It can be shown that carrying 
out the limit in the regularization parameter, the theory in fact becomes trivial in the sense 
that the renormalized coupling constant vanishes. For our purposes, however, it is more 
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adequate to keep the cut-off large but finite, since we want to study A$ 4 -theory only as a 
model for the Higgs. Consequently, we express the quantities A, 4>, m in redefined ones, but 
keep in mind the eventual limit A — ► oo. 

We regroup the corresponding terms fli~4] ) in the following — admittedly suggestive - 
way 

„2 (1 _ = ( _ m 2 + A _A* + A^ 2( _ \x>) + A (— - —(1 + F(k))) 
p v 2 ; v 16vr 2; 2 r v ; \24 8tt v 

and expect that in the renormalized equation, the coupling constant A gets replaced by 
the renormalized one Xr in the rightmost term. Multiplying the whole equation by Xr/X 
suggests the following identifications 

Since we have only calculated to order A we are free to neglect terms 0(X 2 ) which have to 
be fixed by higher order loop calculations. Thus the equation for may be recast in the 
more convenient form 

0(1 - ±V) = 

which reveals an essential fact. From this equation and the redefinition of the coupling 
constant A in fll5D we find immediately that <pX = 4>rXr and consequently k = Kr. The 
invariance of k under the redefinitions is important for the solutions of the gap-equation. 
As opposed to this, if k had a non-trivial scaling with a divergent quantity, the integral 
equation would have only the physically uninteresting solutions for k = or k —>■ oo in the 
limit A — >• oo. Finally we are able to reexpress the gap-equation for the mass parameter in 
terms of the redefined quantities, 

^ = - m % + ~~4>\ + ^ (^ - £(1 + F{k))) . (16) 
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phi 

FIG. 3. The mass parameter fi for the temperatures T c (left) and T (right) for Xr = 0.1 in 
units of rriR. The lowest curve corresponds to the solution \i of the mass equation. The middle 
curve are the values found without the correction term F(k) (superdaisy resummation) and the 
upper curve corresponds to the 'infrared mass'. 

DISCUSSION 



The mass equation (jig) — together with the known function F(k) — may be used to 
determine the effective mass parameter \i for a given value <pR. First we note that the 
contribution F(n) entirely comes from the resummation of the second graph in Fig. 1. 
For small values of the parameter k, the additional term can be neglected with respect 
to 1. On the other hand, for large values of k which in the non-trivial vacuum <Pr ^ 
means small fj,, due to the large k behavior of F(k) ~ k 1 ^ 3 , the mass parameter cancels 



out in the next to leading order contribution in equation fllBD and we are left with a term 
^lTF(k) oc T 4//3 (Ar0k) 2//3 . We emphasize that this term is an important contribution near 



the critical temperature T c = uir^A/Xr because it takes on a non-vanishing value for 
li — > 0. In particular, expressed in orders of Xr, since (f)R,min = 0{1) one finds XrhT c F(k) ~ 



XrT c 4/3 (X 



R9R 



,2/3 



O(Xr) which is of the same order as the term Xr<J) r /2 and thus crucial 



at the critical temperature where the leading term —m R + XrT^/2A in ( |16D vanish. 

Due to this behavior the mass-equation does not necessarily have a real solution for all 
values <pR. In particular, at T c no solutions exists for <pR < rriR\/3F^ A / '(87T 5 ) 1 / 4 = 0.245 . . .. 
For temperatures between T c and T = T c (1 — XrF^ 2 / \Z24ir 5 )~ 1 } the solution consists of two 
disconnected branches which join for T > T to give a solution for all values of 4>r (Fig. 3). 
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The mass parameter may be interpreted physically as a kind of 'ultraviolet' mass appear- 
ing in the propagator (|3]) since the self-energy vanishes for large momenta. As opposed to 
this, in the infrared region k — > the self-energy does have the non-vanishing limit fi 2 A 2 /a 2 
which contributes to the mass in the propagator. The corresponding effective masses are 
shown in Fig. 3. 



CONCLUSION 

We have solved the gap-equation consistently to lowest order for a spontaneously broken 
theory. It turns out, that the Hartree-Fock approximation corresponding to superdaisy 
resummation fails to be consistent. Instead, for small effective mass, we encountered an 
additional contribution in the gap-equation which behaves like 2//3 which cannot be found 
in a perturbative calculation. This contribution of course alters the effective propagator and 
consequently the corresponding effective potential. Since it is exactly the region of small 
mass which is of crucial importance for the behavior of a theory near the phase transition, 
we expect significant changes in the transition of field configurations from the false to the 
trivial vacuum and vice versa. We plan to clarify this point in a forthcoming investigation. 
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APPENDIX: 



The conditions imposed to determine a, A are that equation (|8|) is fulfilled exactly for 
11(0) and II"(0). From the first condition, we find, putting the ansatz fllCf ) into (|8|) 

A2 } j_. 1 nr ,2 , 1 A 2 + ^Kp 



, k dz dp^-, N = p 2 + 1 + a 2 . (Al) 

a 2 _i 7 A^ 2 y a 2 + ^ K ^ + p2 ^ > 

The functions a, A can be approximated consistently for large values of k by a simple power 
law, A(k) ~ A ft 2 / 3 , a(/t) ~ ao ft 1 / 3 . To see this explicitly, we plug the assumed behavior into 
equation ( |Al| ) and rescale the momentum by p = ft 1//3 x which transforms the denominator 
AT into 



N = 1 + n 2/3 L, L 



( A 2 -I- T— 

x 2 + ' ' 2 



x 2 



For large values of k we neglect 1 with respect to the term containing L. Thus the integral 
becomes proportional to ft -1 / 3 which is consistent with the power-behavior of A 2 /a 2 and 



A 2 7 T 2 



( [All ) takes the form 

'i =2 S ix h- (A2) 

u o 

Similarly, from the second condition for il"(0), one finds 

Al 7 ( x 2 (3 + 2L 2 + L 3 ) 2x\2 + L 2 ) 2 \ 

4~{ { 3L^ J' (A3) 

2A t (4A 2 + vr 2 x) 8^(-2a 2 A 2 + 6A 2 x 2 + x 3 vr 2 ) 2 



(2a 2 1 A 2 , + a 2 1 7r 2 2 + 2A 2 1 2 2 ) 2 ' d (2^ + agv^x + 2A 2 ,x 2 ) 3 ' 
The set of equations (|A2| , [A3|) can be solved numerically for the constants Ao, ao, and we find 
a = 4.1140123 . . . , A = 4.8610750 .... 

The function F(n) defined in ([13]) may be calculated analogously by the same reasoning as 
above to have the asymptotic form F(k) ~ Fq ft 1//3 where the constant Fq = 0.996322285 . . . 
is given by 

- , , L-x 2 
F = — I dx- 



7T JO 
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For vanishing parameter k, by the bounds given in (|9]), we have H m ax ~ k and conse- 
quently 

ill k) = k— arctan — 

Reexpressed in the functions a, A this means A(k) ~ ^/Gttk 1 ^ 2 , cl(k) ~ a/12 and ~ 
ttk/6. 

Further values may be read off from the following table. 



K 


io- 2 


IO" 1 1 10 10 2 10 3 10 4 10 5 10 6 


a 


3.491 


3.709 4.995 9.098 19.24 41.58 89.56 192.9 415.6 


A 


0.4358 


1.421 5.210 22.31 104.9 493.0 2291 10636 49370 


F 


0.005397 0.04984 0.3287 1.293 3.701 9.043 20.59 45.49 99.15 
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